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Geometric Polarimetry — Part I: Spinors and 

Wave States 

David Bebbington, Laura Carrea, and Ernst Krogager, Member, IEEE 

Abstract 

A new approach to polarization algebra is introduced. It exploits the geometric properties of spinors 
in order to represent wave states consistently in arbitrary directions in three dimensional space. In this 
first expository paper of an intended series the basic derivation of the spinorial wave state is seen to 
be geometrically related to the electromagnetic field tensor in the spatio-temporal Fourier domain. 
Extracting the polarization state from the electromagnetic field requires the introduction of a new 
element, which enters linearly into the defining relation. We call this element the phase flag and it is 
this that keeps track of the polarization reference when the coordinate system is changed and provides 
a phase origin for both wave components. In this way we are able to identify the sphere of three 
dimensional unit wave vectors with the Poincare sphere. 

Index Terms 

state of polarization, geometry, covariant and contravariant spinors and tensors, bivectors, phase 
flag, Poincare sphere. 

I. Introduction 

The development of applications in radar polarimetry has been vigorous over the past decade, 
and is anticipated to continue with ambitious new spaceborne remote sensing missions (for 
example TerraSAR-X [1] and TanDEM-X [2]). As technical capabilities increase, and new 
application areas open up, innovations in data analysis often result. Because polarization data are 
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spatial and vectorial as well as complex, polarimetry possesses a rich mathematical background 
that includes not only standard electromagnetics but also linear algebra and the algebras of 
Hermitian forms and Lie groups [3]. 

Whilst the aforementioned areas are well established mathematically, the development of the 
formal theory of polarimetry historically has been somewhat unusual, in the sense that it has 
seemed to require some uncommon concepts, such as the pseudo-eigenvectors of Huynen [4], and 
more recently consimilarity transformations [5], [6], [7] have been invoked. There has long been a 
sense amongst some researchers that the theory of polarimetry requires reforming, and a number 
of efforts to do so have occurred over time [8], [9], [10], [11], [12], [13], [6], usually evoking 
controversial responses [14], [15], [16], [17]. What is quite interesting about these episodes is 
that expert opinion is often divided over what precisely is wrong with any new approach. 

Clearly any new attempt along this road must be taken not only with care, but with a view 
to justifying the need for and applications of the new approach. In contrast to earlier suggested 
changes to the formalism which addressed quite small parts of the picture (e.g. Graves [8] and 
Hubbert [12]), we have taken a rather radical view. We believe that it is necessary to start ab 
initio — essentially from the Maxwell equations. We are aware that some researchers in the 
discipline see no essential problem, and hence no need for a reform — and they may feel quite 
justified in such an opinion if they continue to obtain results which work in practice. In response 
to this type of objection, we would counter that it is quite possible that in any given application 
area, a certain mathematical formulation works without any contradiction appearing. Problems 
can occur, however, if one tries to extend a representation to another application area where 
the analysis no longer work; as time goes on, it is a virtual certainty that such attempts will 
be made. The underlying causes of problems may be seen to lie in the very pervasiveness of 
polarimetry in scientific applications. The problem lies in the fact that many application areas 
in polarimetry either involve special cases, or avoid certain special cases — for instance, the 
conventions used in the optical community for ellipsometry assume a fixed plane for reference^, 
so that the transition from what radar polarimetrists would call bistatic to monostatic involve 
no ambiguity. In the radar community, on the contrary, a special convention is normally used 



'since ellipsometry is a specular optical technique (the angle of incidence equals the angle of reflection), the convention is 
to use as a plane for reference the plane spanned by the incident and the reflected beams called plane of incidence [18]. 
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such that monostatic scattering in reciprocal media involves symmetric scattering matrices. But 
in moving to the bistatic case, which is becoming an increasingly topical interest, there may be 
inequivalent mathematical representations, not all of which are equally appropriate: we address 
this very briefly below. 

We are aware that the general polarimetric community is not primarily interested in math- 
ematical exposition per se, and there has been some debate as to whether this is the correct 
forum in which to introduce our ideas. However, in our view, so much is wrong with the current 
state of polarimetric theory that it is necessary to explain first what is the matter with it, and 
preferably to provide an answer. Although this is necessarily rather theoretical, we strive to keep 
our development as approachable as possible. So before moving on to introduce the substance 
of our new approach, we offer a selection of concrete examples that demonstrate that the current 
conceptual framework in polarimetry is lacking, and a cause for concern for the rigor of the 
discipline. It should also be reasonably clear that many of the problems are linked, and must 
relate to something fairly fundamental. We hope that it may be seen as worthwhile to invest 
effort in establishing a workable theoretical basis for polarimetry. 

We note, firstly, a well known identity concerning unitary matrices, namely, 

= U-^ or (1) 

= U-\ (2) 

where U denotes the complex conjugation of the unitary matrix U, and [/^ the matrix transpose. 
This relation in either form appears to have had, and fulfilled, its potential to cause confusion 
in the development and application of radar polarimetry formalism. Let us take three types of 
transformation law that occur in application to matrices in the literature: 

A ^ T-^AT (3) 

A T^AT (4) 

A f-^AT. (5) 

Here, we are not necessarily assuming any particular properties for transformation matrix T. 
The first of these ^ is known mathematically as a similarity, the second dH) as a congruence, 
and the third ([5]) as a consimilarity or conjugate similarity [5]. To understand why each of these 
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might be applied in a given situation, it is helpful first to see what kind of object gets mapped 
to what. 

In case ([3]), if the matrix A operates on a vector x, geometrically a point in a space, to produce 
another vector or point x' 

x' = y4x, (6) 

then the similarity transformation T provides a new representation x of the point x, 

X = Tx (7) 

such that transformed operations on transformed vectors are equivalent to operations performed 
on old vectors. In fact, using ^ and (|7]) 

Tx' = ATx x' = T-^ATx = Ax (8) 

shows that the operation A on x follows the same equation of A on x if A is transformed with 
a similarity transformation. Moreover, the similarity applies equally correctly to a sequence of 
transformations A1A2... 

The geometrical interpretation of a congruence (HI) is rather different. In this case the matrix A 
describes a different type of operation altogether that maps a vector x, geometrically a point, into 
an object u known as covector (see [19]). Geometrically covectors, according to the dimension of 
the space being modelled are lines in two dimensions, planes in three dimensions or hyperplanes 
in general dimensions, defined by N points. In geometry, this operation is called correlation 
[19]: 

u' = Ax. (9) 

A transformation T which provides a new representation x of the point x 

X = Tx (10) 

acts in a different way when applied to the covector u (see [19] for a full description): 

u = T-^^u. (11) 

It is very clear, also from the transformation property that vectors and covectors are completely 
different objects, belonging to different spaces. The congruential transformation, obtained using 
and (dH) 

T-^'^n = ATx, =^ u = T^ATx = Ax (12) 
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is different from a similarity because the thing mapped and the thing mapped to are different 
types of objects. An important example of a correlation matrix A is a quadratic form, that 
describes a quadratic surface, such as an ellipse in two dimensions or an ellipsoid in three. A 
congruence then ensures that in the transformed space, points lying on the original geometrical 
figure will lie on the transformed one. All this appears far abstracted from polarimetry, but 
will ultimately figure prominently in our development of scattering matrices, which lies out of 
the scope of the present paper. However, the form of © is instantly recognizable as that for 
basis transformation of a backscatter matrix (for example [9]) — the only problem is to explain 
why theoretically such a form may be required. In the standard formalism the answer is totally 
obscure. Note, also, that a sequence of correlations AiA2A^.. does not make mathematical sense, 
since the domain and range of the transformations are different types of object. The domain is 
constituted by proper vectors and the range by covectors. 

The third example © which relates to consimilarity is an example of an antilinear transfor- 
mation, in which vectors or points are mapped to a new space of vectors or points that involve 
a complex conjugation to correspond. This type of transformation is well described in the text 
by Horn and Johnson [5] which is regularly cited in relation to this type of transformation. It 
is unfortunate that, as far as we are aware, there has been no fundamental physical justification 
for introducing this into polarimetry. Because of (O, it is obvious that in the case of a unitary 
transformation, consimilarity and congruence are identical. Perhaps because of the prevalence of 
Hermitian forms and of unitary transformations in relation to basis change, the apparatus of quan- 
tum mechanics has had a strong influence over the mathematical development of polarimetry. But 
whilst mainstream quantum mechanics fundamentally involves only unitary operations on Hilbert 
spaces, polarimetry in the sense that it is based on classical electromagnetics, definitely does 
involve non-unitary processes, because scatterers and media in general can absorb. We reiterate 
the point that consimilarity is not a principle. Also germane to the previous point is the concept 
of wave-reversal. It is of course well known that, when using the same coordinate system to 
compare two waves that trace the same polarization ellipse, but propagate in opposite directions, 
the rotation of the complex vector reverses, and that in the analytic signal representation, this can 
be expressed by conjugation of the phasor [6]. Consequently, in the Jones vector representation 
[20], the vector components are conjugated for wave reversal. This is sometimes mistaken as 
a justification for adopting the consimilarity formalism in radar backscattering. However, it is 
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quite clear that backscatter is not an antilinear transformation, since no time-reversal occurs. 
Mathematically, the operation of conjugation cannot be used selectively — a point that can be 
missed by the unwary when analytic signal representation is used. Specifically, for the analytic 
representation of a harmonic plane wave, we can write 



Qj(ujt~kz) = QjLUt Q-jkz = Qj{kz~u)t)^ 



where the overbar Qi{^t-kz) denotes the complex conjugate operation, uj is angular frequency and 
k is the wavenumber. The final conjugated result has the same phase velocity as the unconjugated 
wave. Mathematical conjugation cannot be invoked to reverse just space or just time. 

One of the most significant and useful developments in polarimetric analysis has been that of 
the target vector covariance theory introduced by Cloude [21], and which has been developed 
to a mature state [22]. Here, the starting point is a linear representation of scattering matrices 
S using the Pauli matrix decomposition: 

S = ao-Q + bai + ca2 + dcr^ = (14) 

+ c 





where ai are one matrix representation of the Pauli spin operators. In this, the representation 
of the scattering matrix is taken as a given, and the group-theoretic techniques that have arisen 
in this area arguably have more to do with information theoretic principles than the physics of 
scattering. If, however, one wishes to connect this theory to the underlying physical principles, 
problems can arise. It was noted in [23] that the vectorial interpretation of scattering matrices can 
be better related to the Poincare sphere if the original Pauli matrix decomposition is modified. 
This problem can be related to the distinction between Q and (21) . Mathematically, Pauli matrices 
ai are defined via multiplicative properties 

Cri(rm = SlmC^O + j ^ImnCTn, (/, m, ?2 = 0, 1 , 2, 3) (15) 

where dim is the Kronecker delta symbol and eimn is the antisymmetric Levi Civita symbol 
[24]. The relation (fTSi) means that the actual realizations of the matrices are just one example 
of an equivalence class that transforms according to a similarity transformation. Geometrically, 
applying the congruential transformation & to the Pauli matrices makes no sense, and yet, under 
a basis transformation, (Q) is precisely what is required. It is also a main plank of the target 
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covariance theory that it extends to bistatic scattering by an obvious generalisation, by including 
the antisymmetric Pauli matrix 

Geometrically, however, the distinctively different Pauli matrix is ctq. As [23] shows, an im- 
portant simplification arises if a twisted representation is used — and in fact this can be easily 
related to the concept of Huynen's pseudo eigen-equation [4]. Target covariance theory 'works' 
until one tries to relate it to the geometry of the Poincare sphere. A recent example where a 
misunderstanding of the underlying mathematics has resulted in a very confused picture was 
in [25] which attempted to apply a quaternion analysis to the bistatic case, but neglected the 
same point. It should in fact be obvious that, within the backscatter alignment (BSA) convention, 
multiplication of scattering matrices makes no physical sense (for it cannot be used consistently 
to describe multiple scatter). Modelling scattering matrices on a division ring (the technical 
designation of Pauli matrices [26]) makes no sense without a fundamental modification. 

To summarize the first part of our introduction, we have shown with a few examples that the 
mathematical and conceptual framework that pervades much of the radar polarimetry literature 
rests on very shaky foundations. Even if much of the analysis turns out correctly, both the physical 
and mathematical principles on which it rests have become obscured by a great deal of faulty 
reasoning over the last half century. We now commence the establishing of a new framework, 
which ultimately should unify the treatment of coherent and incoherent polarimetry in traditional 
forms and also naturally extends to the concepts of target covariance. The anticipated growth in 
bistatic applications suggests to us that it is timely to present such a unified theory, particularly, 
since it is capable of handling polarization states for waves propagating in arbitrary directions 
in the three dimensional space. We also hope that, once and for all, the origin of apparent 
idiosyncracies of polarimetric algebra will be settled. 

Before going on to the main exposition, it is appropriate, given the novelty in the eyes of the 
polarimetric community of the formalism we employ, to say a few words concerning the scope 
of our objectives here and later, and the importance of what is gained. The most important single 
aspect of this work is emphasised by the first word in the title, which is geometric. This end is 
far more important than the means, which is the use of spinor algebra. If the ideas could have 
been presented in some other way, we would gladly have done so, however, it turns out that 
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spinors are so perfectly adapted for representing polarimetric objects that we would otherwise 
have had to re-invent them. Spinors were introduced by Cartan [27] in a geometric context, and 
may perhaps be best summed up as doing much the same for geometry as imaginary numbers 
do for algebra. It may be remarked that spinors appear rather similar to Jones vectors, and the 
question may well arise, as to whether we are merely re-inventing them. Although these objects 
appear virtually identical in some contexts, there is a problem with Jones vectors in the sense that 
they seem to have ambiguous properties. In one sense they are extensions of Euclidean vectors, 
although they also may be subjected to unitary (basis) transformations. In orthogonality relations, 
Hermitian products are employed, whilst in antenna voltage equations a Euclidean inner product 
is appropriate. Mathematically, this ambiguity makes it difficult to assign any unique property, 
and it follows that the mathematical nature of the polarimetric model as a whole is hard to pin 
down rigorously. What is fundamental in this our first paper is the assignment of a coherent 
polarimetric state to a spinor, as a geometric entity, which turns out to be a line in projective 
space. Conjugate wave states are represented by conjugate spinors, and their intersections define 
Stokes vectors. Because spinors are truly geometric entities derived from the electromagnetic 
field, we are led to a quite outstanding result, that we are able to identify the Poincare sphere with 
the sphere of normalised wave-vectors in the Fourier domain. Thus, no longer are polarization 
states (in the form of Stokes vectors) consigned to some abstract space. Rather, we can construct 
objects all in the same space. Given this, we can represent arbitrary geometric rotations in the 
same analytical framework as other polarimetric relations. 

The key to achieving this will be the introduction of a hitherto 'hidden' object, an implicit 
reference spinor, which we name the phase flag. This step should enable physicists to overcome 
their natural reticence to describe polarization vectors as spinors, which are in this algebraic 
form associated in field theory with 'half-spin' particles like neutrinos. Physicists' explanations 
for the SU(2) representation of light polarization (e.g. [28]), tend to appear half hearted and 
lacking conviction. In essence, when a reference spinor is taken into account, the extra half-spin 
required by a photon is accounted for. The amazing simplicity of this idea seems not to have 
surfaced in the physics literature, which probably accounts for why spinors have to date figured 
little in polarimetry in any formal sense. This also explains why Stokes vectors do not transform 
relativistically [29], even though they appear to be rather like four vectors. 

Our presentation is rigorously derived from Maxwell theory for harmonic waves, and by the 
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use of a projective formalism that extends to the spinor case [30] provides a complete and 
consistent geometric theory of polarization. The importance of a geometric interpretation for all 
the objects in polarimetry can hardly be overstated. For once one has this, algebra is only an 
adjunct; geometric relations remain true whatever basis transformations are applied. 

This paper emerges after having worked out the implications of our ideas to a wide range 
of polarimetric concepts. Owing to the space required to set out the fundamentals from first 
principles, we limit ourselves here to the modest objective of introducing the derivation of the 
formalism as it applies to polarization states. We have to defer its application to scattering 
matrices and other entities which we hope to address in later papers. We can however assert that 
the geometrization of polarization leads to very simple and elegant statements about scattering 
matrices and their properties in a most convincing way. Ultimately, it will be seen that much of 
the algebra appears quite similar to that in conventional polarimetry, but the key point is that 
there is in the end a clarity of concept and consistency in algebraic characterization of objects 
that does not currently obtain. 

II. Spinors and geometry 

The concept of polarization state has long been associated with the Poincare sphere in a 
construction sometimes known as a Hopf mapping through the polarization ratio [31], that is, 
fundamentally, a geometric construction. The Hopf map [32] takes the points (a, b, c, d) on the 
unit sphere in M^, labelled as = {{a,b,c,d) : + b"^ + + = 1} to points on the unit 
sphere in M^, labelled as S'^ = {{x,y, z) : x"^ + y"^ + z"^ = 1} through the following relations: 

x^2{ac + bd), y^2{bc-ad), z ^ + b'^ - - d'^. (17) 

Introducing a — a + jb and P — c + jd we have 

+ b^ + + d^ = 1 ^ |ap + |/3p = l. (18) 

We can recognize three of the Stokes parameters in {x,y,z) and the components of the electric 
field in (a,/?). Since 

f = x + iy = ^ = i±^ = ^. (19) 

1 + z a + jb a 
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in polarimetry the Poincare sphere is recovered through the stereographic projection, where the 
polarization ratio, C, = is first mapped to a point on the real Argand planer and then the 
mapping is normally performed through the stereographic projection from the Argand plane to 
the unit sphere^ as shown in Fig. [TJ The relations (l20l) become: 

C + C j(C-C) 1-CC .... 
X = — , V = — = , z = — . (20) 

CC + i CC + i CC + 1 

and they codify the Poincare sphere. 

The aim of this paper is to perform a construction of the Poincare sphere as a collection of 
complex lines and not as a Hopf map. Each coherent state of polarization will be one of the 
complex line generating the sphere, obtained in form of a spinor from the full electromagnetic 
field. Such a line can be parameterized by a complex number equivalent to the polarization ratio. 
The first step is now to see how a spinor can be interpreted as a generating line of a sphere, 
and to see the links between a sphere, a complex line and a spinor. We will attach the physical 
meaning to the spinor in the next sections to establish that this sphere is the Poincare sphere. 

A suitable starting point is the phase of the wave. The analytical representation of a plane 
wave is proportional to [34] 

giM-kx) (21) 

with uj the angular frenquency and k = {kx,ky,kz) the wavevector. According to relativity theory, 
the phase ip = tut — k ■ x is invariant [34]. This means that in two different reference frames 
(two observers in uniform relative motion) the plane wave would have different frequency cu 
and wavevector k but the phase Lf would be the same. As a consequence, the invariance of the 
phase corresponds to the invariance of a sort of scalar product between two vectors with four 
components —k^, —ky, —k^) and {ct, x,y, z): 



(^^,-k^,-ky,-k,Y {ct,x,y,z), (22) 



^where the real part of is represented by a displacement along the X— axis and its imaginary part as a displacement along 
the Y— axis of the Argand plane. 

^ The mapping can be performed in any polarization representation, since it is a function of a complex variable. In this case 
and throughout the paper, we use circular basis, following the notation of Penrose [33]. The polarization ratio is chosen as 

— where the North pole has a polarization ratio of infinity (a = 0) representing left-handed circular polarization and the 
South pole has a zero polarization ratio (/3 — 0) representing the right-handed polarization. For a good summary of polarization 
ratio in different polarization basis see Mott [31]. 
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where c is the speed of light in a vacuum. Because of this invariance, the frequency and the 

wavevector of any plane wave must form a 4— vector. The tensor description elegantly expresses 

the invariance property and highlights the different behavior of the two 4— vectors in case of 

change of reference frame. In fact, the 4— vector {ct,x,y,z) = is a true vector since its 

dimension is a length, instead the wavevector is a 4— gradient of a scalar invariant 

/cu, , ,\ f I d d d d \ ^ 
-, -K, -ky, -k,) = [ - — , ] V = Oa^ (23) 

V c / \cotoxoyozJ 

n 

which has the dimension of the inverse of a length. The 4— vector x'^ transforms contravariantljO 
with respect to the 4— vector ka = —kx, —ky, —k^). The simplest demonstration of this is that 
if we change our unit of length from meters to centimeters, the numerical values of x"" would 
scale up while those of ka would scale down. The tensor notation automatically takes account 
of this and tensor jf| like x"' are called contravariant while tensors like ka are called covariant. 
The tensorial notation of the product (|22l) would be: 



^ ^ kaX — kaX , (24) 



a=0 



where we have dropped summation sign adopting the Einstein summation convention where 
upper indices {x"") are paired to lower (A;^)^. 

We shall try now to obtain representations for the complex lines generating the sphere. In order 
to do so, we consider now the contravariant 4— vector = (r = ct,x,y,z). It is isomorphicj^ 
with a Hermitian matrijqj which may be parameterized in the form 



r, x,y,z e M. (25) 



T + z X + jy 
X — jy T — z 

The condition that the matrix be singular may be expressed as 

- x'^ -y^ - z'^ = 0. (26) 

''which means that the inverse transformation has to be used for ka. 
4— vector is a type of tensor. 

^Inversely, if the distinction between contravariant and covariant were not made, scalar products would only be invariant under 
isometric transformation such as rotations or reflections, and not under unitary transformation in general. 

'An isomorphism (from Greek: isos "equal", and morphe "shape") is a bijective map between two sets of elements, 
matrix M is Hermitian if M = where denotes the conjugate transpose. 
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In special relativity this condition is satisfied by points lying on a light-cone through the origin 
[35]. An alternative interpretation is for (T,x,y,z) to be considered as homogeneous projective 
coordinates [19]. Let us suppose we have a point (X,Y,Z) in the Euclidean space M!^. To 
represent this same point in the projective space P'^, we simply add a fourth initial coordinate: 
(1, X,Y, Z). These coordinates are the homogenous coordinates of a point in the projective 



space P^. Overal 

(a, aX, aY, aZ 



scaling is unimportant, so the point {1,X,Y,Z) is the same as the point 
In other words, 

(r, X, y, z) = (ar, ax, ay, az). (27) 



Since scaling is unimportant, if the coordinates (r, x, y, z) are considered the homogeneous 
coordinates of a point in the three dimensional projective space P'^, the equation (|26|) then 
defines projectively a sphere (or more generally, a quadric surface): 



x\2 fy\'^ fz^'^ 



This is reminiscent of the reduction of polarization states of arbitrary amplitude to a unit Poincare 
sphere. 

From one side the vanishing of the determinant of the Hermitian matrix (l25l) allows us to 
build projectively a sphere but also, since it is Hermitian, to express it [24] as the Kronecker 
product 

I T A- 7 nr A- in \ 

(29) 

A'e' ee') 

where (g) denotes the Kronecker product, and the overbar denotes complex conjugation^. The 
matrix on the second line in (|29l ) built with the complex 2— vector (^°,^^) is in fact singular 
and Hermitian. The basic ideas here will appear familiar in polarimetric terms as relevant to the 
construction of the coherency matrix [36]. 




'For any a ^ 0, thus the point (0, 0, 0, 0) is disallowed. 

'"The reason for labelling the conjugated elements with a primed index will become clear in the sequel section IIII-BI 
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We can notice that the condition necessary to build a sphere projectively is the singularity of 
the matrix and not the Hermiticity. In fact, if we relax the Hermitian constraint we have that in 
general {T,x,y,z) are not real, but they still build a singular matrix and projectively a sphere. 
The points {T,x,y, z) are complex points lying on the 'real' sphere defined as 

-x"^ -y^ - z^ = 0, r, x,y,ze C. (30) 

Since the matrix is not Hermitian anymore, we have to redefine the terms of the Kronecker 
product (l29l) 

/ T -i- 7. T -I- 17/ \ 

(31) 





For general complex pairs (^°,^^), {rf.rf), the matrix is still singular but not Hermitian. Now 
if {rp , rf) be fixed, it can be seen that points 

{T,x,y,z) (32) 

form a linear one dimensional, complex projective subspace, (^°,(^^) being variable. We want 
to show now that this one dimensional projective space is a complex line generating the sphere. 
This is the main result of this section. We will use this conclusion to define the polarization state 
and the Poincare sphere in the next sections. It should be clear that the same argument may be 
applied keeping (^°,(^^) constant and varying [rf.r]^). In this way, a second line is generated. 
Every single independent (^°,(^^) {rp^ff) pair results in a unique point of the sphere. 

The singularity of the matrix (|3TI) implies the existence of a sphere in the projective space 
We want now to see what is projectively the complex 2— vector {^^,^^). As usual, we consider 
(^°, as homogeneous coordinate, which means that (1, C = or any non-zero multiples of 
it, represents projectively the same vector. So we associate the vectors ((^°, (^° 7^ 0) with all 
the complex numbers together with the infinity (0,1) and we obtain a space which is usually 
labelled as CP^ [19]. Regarding the dimension of this space, this depends on the point of view. 
From the perspective of the real numbers already C is a two-dimensional object often called 
the Argand plane0- As we have already noticed in polarimetry theory this is the traditional way 



"The space C differs from CP^ by one point, tfie infinity. 
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to arrive at the Poincare sphere, where C might be interpreted as a polarization ratio. However, 
there is the other perspective, the complex perspective, from which C is just a one dimensional 
space. It contains just one complex parameter. In the same way as M is the real numbers line, 
one would consider C as the complex numbers line and so CP^ is called the complex projective 
line. The elements of CP^ are the complex 2— vectors (^°, ^^), known as spinors and represented 
by a complex line. The point P on the line is represented by the projective parameter C = |o 
once two reference points R and Q have been chosen (see Fig. O: 

p = ^Oi? + ^iQ, ^ P = R + (:Q. (33) 

At this point we have a complex projective line and a projective sphere linked by the relation 
(|3T| ). We want now to discover the geometric significance of this link: it illustrates an example in 
the theory of projective geometry that through any point on a quadric surface there pass two lines, 
each of which lies entirely within the surface [37]. To those unfamiliar with complex geometry it 
may appear surprising that a sphere contains straight lines: this is a fact that escapes us because 
on a sphere or any quadric with positive curvature only one point on each such line is real; all 
other points are complex. A simple way to "see" how a complex line can belong to the sphere 
is to consider planes intersecting the sphere. We consider for simplicity a sequence of planes 
parallel to the x^y— plane, cutting the 2— axis at Zf.. The equation in the plane of intersection 

x^ + y^ = l- zl (34) 

is the equation of a circle which degenerates in a point of tangency {zc = ±1). In this case, the 
point of tangency is not the whole solution. In fact, we have also two intersecting complex lines 
with gradient ±j 

+ y2 = ^ X = ±jy (35) 

having one real point (1, 0, 0, 1) at their intersection. What is surprising is that the two lines lie 
completely on the sphere surface. In fact, each point on these lines has coordinate (1, ±jy, y, 1) 
which belongs to the sphere. This is a simple example to illustrate that a sphere contains complex 
points (|30l) and in particular lines built up with complex points ! Now, the next step is to see how 
such lines can generate the sphere surface. Since it is difficult to imagine this we can consider 
instead a quadric surface with negative curvature^. For such quadric surfaces, the generating 



'"In reality, in tiie complex projective space all the non-degenerate quadrics are indistinguishable from each other [19]. 
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lines can be wholly real, a fact that is exploited architecturally, e.g. in the design of cooling towers 
as cylindrical hyperboloids [38]. In Fig. [3] we show how a line can generate an hyperboloid. 
We can easily see that i) any point of the generating line is a point of the surface, ii) there are 
two families of generators, and through each point of the surface there pass two generators, one 
of each family, iii) generators of the same family do not intersect, iv) each line of one family 
intersects with every other line of the other family. 

The point of the preceding discursive outline is to emphasize the central principle of geometric 
polarimetry, which is to identify the complex 2— vectors such as 

^^ = ie,e) (36) 

with one of the two sets of complex lines generating the sphere. The complex 2— vector is 
known as a spino^ [39], [27], [40]. We should stress at this stage that the sphere in question 
is to be thought of as the sphere of real unit vectors in three-dimensional space. The concept of 
Poincare sphere will be derived later. 

In order to establish the notation that we will use in the sequel of the paper, since the singular 
matrix (|3TI) is the outer product of two spinors and the conjugate of r]^, we can label it as 
follows: 



^00' j^oi' \ I r + z x + jy 



=X^'' = I 1 = 1 "'^ I • (37) 

We can notice that X"^^' is considered as a spinor with two indices; generally spinors can, like 
tensors have any number of indices. 

We have started from the phase of a plane wave built up with two 4— vectors, x"" and ka- But 
so far we have considered only the contravariant coordinate 4— vector in the physical space. 
Now since ka is also a 4— vector we can think to make the same considerations we have done 
for x"". However, there is a difference. The 4— vector ka is covariant. Contravariant and covariant 
4— vectors are called duals of each other. We want to explore the meaning of duality in the 
projective context. First, we consider the equation of a plane in the three dimensional Euclidean 
space: 

/(x, y, z) = ux + vy + wz + g = 0. (38) 



'^^The term was originated by Cartan [27]. 
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The normal vector to the plane is given by 

N = V/ = {u,v,w). 



(39) 



The plane may be characterized via its normal vector from the origin, so a set of three coordinates 
{u,v,w) can refer to a plane rather than a point, giving rise to a 'dual' interpretation. If we 
consider the projective space, we add a fourth initial coordinate and we have that the condition: 





T 




u 




X 


V 




y 


[w) 




v) 



(TT + ux + vy + wz = 



(40) 



may be seen as the condition for a plane {(7,u,v,w) to pass through the point (T,x,y,z). 
Comparing (l40l) with (|38l) we have that q = ar. Since the set of homogeneous coordinates 
{a, u, V, w) are defined as a gradient (|39l ) we associate it with a covariant 4— vector and we 
associate the set of coordinate (r, x, y, z) interpreted as a point with a contravariant 4— vector 
x". The linear equation of a plane (l40l) can be written as g^x" = 0, with qa = (a, u, v, w). Hence 
the components of any covariant vector qa are to be regarded as the coordinate of a plane. 
Reconsidering the phase of the wave (l23l) . the zero phase 

/ ^ V (r\ 

C 







(41) 



X 

-ky y 

in the projective representation states that the plane ka passes through the point x'^ which belong 
to the sphere 1 - (f)^ - (f)^ - (f)^ = 0. 

At this point, we can make the same consideration we have done for x"" keeping in mind that 

ka = (7, -K, -ky, -k-J) is a plane. 



The singularity of the Hermitian matrix 

Kab' = 

defines projectively a sphere 

kx 
ko 



IKoo' 






Uio' 







ko kz knc + jky 
'kx ~~ jky k() + kz 



(42) 



(43) 
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with ko = ^. The sphere this time is an envelope of tangent planes ka- The singularity of the 
matrix (|42l) can be expressed as: 



kan^'^kb = 





A 








o\ 




ki 





-1 








h 


k2 








-1 





k2 












-1/ 





(44) 



where ki = —k^, ^2 = —ky, k^ = —kz, so that from here on numerical as well as symbolic 
index positions are to be explicitly interpreted as contravariant x"- or covariant ka according 
to the position. We shall refer to ^l"^ as the wave sphere, since (|44|) expresses the Fourier 
transform of the wave equation in free spaceLJ. Now we can also give a projective meaning 
to k"- = Qf^^kb = (fco, k^ = kx, k"^ = ky, k^ = k^) as the point of tangency of the plane k^ on 
the sphere Qf"^. In the Fourier domain, the interpretation of the 4— vectors ka as homogeneous 
coordinates in three-dimensional projective space implies normalization of the frequency 

V kQ ko koj 

As we consider harmonic waves individually, and in many applications a quasi monochromatic 
assumption is justified this creates few problems. In return, the interpretation of the linear algebra 
as three-, rather then four-dimensional is beneficial from the point of view of visualization. 

As for X^^' in (|37] ), the singularity of the matrix (|42l) . allow us to express it as the Kronecker 
product 

Kqi' \ I Kq 

KlV I \Ki 



Kab' 




(46) 



this time with the covariant spinors 

k,a = {ko,Ki) and Aa=(Ao,Ai). 



(47) 



Like the contravariant spinor (|36l ) they also represent complex lines on the sphere. The concept 
of duality is still valid for spinors, namely covariant and contravariant spinors are still duals of 
each otheo. However, in the three dimensional projective space P'^, the dual of a line is a line 

'"*As we will see in the section ITV-B I in the Fourier domain da jka- 
'^It means that they transform covariantly to one another (see footnote |4j. 
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or lines are self-dual. In fact, a line can be built linking two points but also it is the intersection 
of two planes, duals of points. Given two points p"" = (p'^ , , p'^ , p^) and q'^ = {q^, g^, g^, g^) in 
homogeneous coordinates, the projective description of the line passing through the two points 
are given by the numbers: 

p'q'' 



p'q^ 



which build the tensor 



lab 



p^q^ - q^p\ 



(48) 



(49) 



Since / = — I the tensor is clearly skew-symmetric: 

/ 
/ 

and therefore the distinct elements are reduced to six 
However, they are not independent since they always satisfy 

/10/23^;20^31^ ^30^12 _Q 








_;20 


-/30\ 


10 









20 









30 




-/23 


/ 



(50) 



(51) 



(52) 



which is the determinant of a 4 x 4 matrix [jf , q°- ,p°- ,q°-) that is identically zero. The coordinates 
(|5TI) connected by the relation (|52|) are called Pluecker (or Grassmann) coordinates of a line. 
Again overall scaling is unimportant, namely the set 



(20 



(30 



23 



f31 



(53) 

represents the same line as (|5T1) does. If the first coordinate of the points is not zero, it is easy 
to show that the coordinates have a nice Euclidean interpretation, namely 

(/10,/20^30N 

^ ^ (54) 

(/23,/31^;12)^pXq 

12 3 12 3 

with p = ^, ^), q = ^) and where x denotes the cross product. The first set of 

coordinates describes the direction of the line from q to p and the second describes the plane 
containing the line and the origin. The condition (l52l) is equivalent to the identically null product 

(p - q) • (p X q) (55) 
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(56) 



where ■ denotes the scalar product. 

Alternatively, we consider the planes pa = VtabV^ and qa = ^ab'f- We define the skew 
symmetric tensor lab 

Lb = PaQb - QaPb 

whose components are related to the components of /"^ in ([49] 

/ /2o 

-/lO 
y_/30 ^31 _^23 



ab 



/30 \ 

/ 



(57) 



which is the dual of the tensor *rab representing the line intersection of the planes pa and Qa'. 

*rab = Lb. (58) 

The dual of a tensor is defined through the full antisymmetric Levi Civita symbol e'^^^'^ and it 
is related to the tensor r"^ : 



, ab 



1 

2 



(59) 



which components can be written in function of the components of 



, ab 












23 





_;30 




31 









12 


-/20 




/ 



(60) 



(61) 



The Pluecker coordinates of the line intersection of pa and qa will be the set: 

{^23, ^31, ^12, —^10, ~ho, — ^3o}- 

Again considering the Euclidean interpretation as in (|62l) . the vector p and q now represents the 
normal to the planes. For this reason 

^ ^ (62) 

(_/10^_/20^_;30) _q_p 

this time the first set of coordinates namely the direction of the line is described by p x q and 
the second set namely the plane containing the line and the origin is described by q — p. We 
emphasize that this is only an Euclidean interpretation that can help us to visualize things but 
the Pluecker coordinates are coordinates in the projective space and not in the Euclidean space. 
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III. The algebra of spinors 

In the previous section we introduced contravariant and covariant spinors by interpreting them 
as representations of the complex lines generating the sphere. In this section we will discuss the 
algebra of spinors especially the relation between covariant and contravariant spinors and the 
conjugate spinors. 

Spinors were first introduced by Cartan [27] in 1913 and later they were adopted in quantum 
mechanics to study the properties of the intrinsic angular momentum of the electron and other 
fermions. Today spinors are used in a wide range of branches of physics and mathematics. 

This paper follows the notations and the conventions of the book by Penrose and Rindleif^ 
[33] which provides a very clear idea of what spinor representation signifies. 

The spinor has four basic forms each of which follows a separate and distinguished 
transformation law. This is the strength of the spinor notation, especially for polarimetry: the 
labelling of a spinor of one of the four type immediately dictates its transformation properties. 

Many of the equations of polarimetry involve nothing more complicated than 2x2 complex 
matrices. By understanding the significance of spinor index types, it becomes clear that these 
are potentially 4x4 different characters of 2 x 2 complex matrices; little wonder therefore that 
confusion may arise over the correct transformation laws! 

We start now to describe the four forms of a spinor and their transformation law, and later 
after the derivation of the polarization state from the electromagnetic field tensor we will show 
why and how the distinction of the same object in four forms is fundamental in polarimetry. 
The four spinor types can be categorized as: 

• the contravariant spinor 

• the covariant spinor ^a, 

• the conjugate contravariant spinor , 

• the conjugate covariant spinor ^a'- 

All these different types have definite and different physical or geometric significance, and must 
not be confused with one another. 

In the next two sections we describe how to raise and lower indexes and the conjugation of 
a spinor. The significance of raising and lowering is to make a relation between one type and 

"which are the same as the conventions of the book of Misner [24]. 
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(63) 



another. This can be seen as a mapping between objects that are dual to one another through 
some fixed geometric property. 

A. Index raising and lowering 

In this section we introduce the metric spinor cab = e"^^ which tells us how to compute 
distances in spinor space. We will see that it is the object responsible for lowering and raising 
the indexes and which can test the independence of two given spinors. 

We have seen in the previous section that projectively the following components of a spinor 
are equivalent: 

'e\ / 1 \ /I 

The spinor components may be interpreted as the homogeneous coordinate of a point P 

on a complex line (see Fig. [2l). Then the complex number ( is the inhomogeneous coordinate 
of the point on the complex line known as the qffine coordinate of a point. Following this 
terminology we can see that the difference between complex numbers ^, ij, namely the affine 
distance between the two points P — P' may be expressed in terms of their homogeneous 
coordinates as 

Now reinterpreting the same relation using the homogeneous coordinates we can recover the 
invariant skew-symmetric spinor eAB called metric spinor. 



(64) 





AB 



(65) 



where we adopt as usual the Einstein summation conventioit_j. The metric spinor has components 



<^AB 



ei2 



-£21 — 1, eii — 622 — 



(66) 



namely it is a completely antisymmetric spinor {eab = —^ba)- 



"The expression eABC^rj^ stands for 5I1a=o J2b=o '^ab^.^'O^ 
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Since the distance is always related to some kind of product, we can symbolically introduce 
the inner product 

eAB^^v"" = CbV^ (67) 



where the covariant spinor 



?B = e^s^^ ^0 I ^^^^ 

is obtained from the contravariant one with the lowering operator gab- Here, as usual, we use 
the Einstein summation conventioi^f . 

Likewise, e^^ which has numerically the same element values as gab is a raising operator, 
performing the inverse mapping from covariant space to contravariant: 

= e^^^B. (69) 

We note that 

eAB e^"" = (70) 

is the Kronecker delta symbol 6^^ in spinor space. 

It may also be recognized that the product (|67l) vanishes if the ratios of the spinor coordinates 
of ^ = 1^ and ^7 = ^ are equal, namely 

CBr/^ = ^ e^'ocr/^. (71) 

Geometrically, if the two points on the line coincide, the distance vanishes. So, we have arrived 
at a symbolic mechanism for testing the linear independence of two spinors: 

UV^ + 0. (72) 

We want finally to show how spinors transform for linear transformations. Since linear trans- 
formations map contravariant spinors to contravariant, we can express the transformation law 
for a spinor as 

e - ~e = i^\e- m 

"*which means = I]a=o '^AsiJ^ , which can be read = cooC" + ^loC^, (,\ = eoi^° + eiiC^- 



23 



with L^^ the spinor describing the linear transformation. Using their matrix representation we 
can write: 



(74) 



where is represented by a 2 x 2 matrix. 



B. Conjugate spaces 

In the last section we have introduced two types of spinors, the contravariant and the co variant 
and their transformation law. The other two kinds of spinors are the 'primed' ones belonging to 
the conjugate space. Thus, 

V = V^' =\ I (75) 




in which the numerical values of the components are conjugated, and the abstract index becomes 
primed. The conjugate spinor fj^' will transforms according to the conjugate law 

r = L^B'n^' (76) 

where L^'^, is the conjugate transformation. The operation of priming applies equally to covariant 
spinors, so we have obtained the four categories of spinors: 77^, rjA, f]^' , fjA'- 

Using the explicit transformation laws for the four types of spinors we can infer the trans- 
formation law for spinors with more indices and see why attaching the prime in the index is 
necessary. 

In section |II] we have considered matrices expressible as a Kronecker product of two spinors 

m 

Now using the transformation law (1731) and (1761) for and for f]^' we can show that X^^' 
transforms not simply as a matrix but 

j^AB' _^ j^AB' ^ ^A~B' ^^A^^C-D'IB'^^^ 

= L^X"^^' L'^'o,. (78) 
Expressing the relation below with the matrix representation we can rewrite it as: 

X X = LXL^, (79) 
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where X is the matrix representing X^^' and is the conjugate transpose of the matrix 
representing L^^. 

We can easily see that the prime attached to the index is very essential in order to keep trace 
of the fact that they transform with the conjugate transformation. Clearly, since conjugation is 
an involutory operation, it means that a double primed label reverts to the unprimed. 

Through an example we can see that if L is a rotation in space, X^^' (and consequently 
the 4— vector) rotates by double the angle of rotation of a spinor. Since the group of unitary 
transformation SU(2) is the universal covering group of the orthogonal rotation group S0(3), this 
means that there is a correspondence between elements of SU(2) and of S0(3)[!^. The unitary 
transformation can be expressed with the Cayley-Klein parameters in terms of the Euler angles 
if), 6 and ip for the rotation considered [41]: 

COS ^ e-' 2 7 sm ^ e-' 2 \ 

, • (80) 

J sm I e 2 cos | e 2 / 

The example is a rotation in space around the z direction by angle a. The unitary spinor is 
in this case 



e-^2 
e^t 



< = r A (81) 



According to ([73]) . the new spinor ^ after the rotation will be: 

i^^.^e-i^- (82, 

For a = 27r this rotation brings back to the original value X^^' (according to (|78] )) but reverses 
the sign of the spinor 

A MATHEMATICA package [42] has been developed in order to perform symbolic calculation 
with spinors and tensors for radar polarimetry using the Penrose-Rindler notation [33]. 

In our view, the correct attribution of covariant, contravariant, and primed or unprimed is 
crucial to a proper understanding of polarimetric transformation properties. A particular source 
of confusion in polarimetry arises from a restricted symmetry, in that for any spinor has 

"The correspondence SU(2) S0(3) is a 2 : 1 one [33]. 
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numerically the same components up to a phase factor as the 'orthogonal' spinor. In fact, two 
independent spinors satisfy (1721) . In particular, if 

UV^ = 1, (83) 

the pair {^"^,r]"^} forms a spin frame which allows any spinors to be represented in component J^. 
The spinor r^^ then will have component^l 

Comparing them with the components of the spinor ^a', obtained using (l68l) and the co variant 
form of (1751) we obtain the components proportional to rj^: 

FO' 



U' = I I . (85) 

This symmetry is equivalent to the property of unitary matrices that their conjugate transposes 
are their inverse. This conflation of properties of different character spinors (which the Jones 
vector calculus fails to distinguish) means that if one restricts consideration to unitary processes 
the statement of the transformation rules in the Jones calculus is not uniquely determined. An 
example of misinterpretation is that a backward propagating elliptical state of polarization is 
considered conjugated. For example this is implicit in Graves paper [8]. Hubbert [12] accepts 
this as part of the backscatter alignment convention. Lueneburg [6] argues its legitimacy by 
the use of the time reversal symmetry, when it holds. However, some radar problems involve 
propagation in lossy media. If backscatter with propagation through absorptive media is to be 
treated within the theory then confusion over this accidental symmetry must be avoided. 

IV. Maxwell's equations and the wave spinor 

In the polarimetric literature the state of polarization is customarily described with reference 
to the electric field vector, although in earlier literature it was the magnetic field vector. Notwith- 
standing, to divorce one from another in an electromagnetic wave is artificial since one can never 
propagate without the other. 

^"A spin frame constitutes what we usually call a basis. 

^'if we assume that has unit Hermitian norm, namely if = 1- 
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We will start from the electromagnetic field tensor which contains both the components of 
the electric and magnetic fields. We will extract from it the direction of propagation to obtain 
the electromagnetic potential. This will be the right quantity to use to derive the polarization 
state since we will be able to establish from it the form of spinor with one index that is 
usually designated as the coherent polarization state without recourse to the complefixication 
of the Euclidean electric field E. In this way the nature of complex unitary transformation for 
rotation in space will become very clear. 

In the next section we derive the spinor form of the electromagnetic field and in the following 
one we project out from the electromagnetic potential a spinor containing the polarization 
information. 



A. Maxwell's equations and the tensor and spinor forms of fields 

Whilst in engineering and applied science, the vector calculus form of Maxwell's equations is 

prevalent, they are more concisely formulated in tensor form. The electromagnetic field tensor 
Fab contains the components {E^, Ey, E^) of the electric field E and {Bj.-, By, B^) of the magnetic 
induction B [34]: 
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(86) 



(87) 



and a further field tensor Gab for linear media containing the components {Dx, Dy, D^) of the 
electric displacement D and {Hx, Hy, Hz) of the magnetic field H, 

/ 



G. 



ab 
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(88) 
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the Maxwell equations in S.I. units can be written as: 

V B 



* TTiab 



daG 



ab 







J" 



VxE=-f 



VD 



P 



V X H = J + 



m 

at 



(89) 



(90) 



where J'' = {cp, J^, Jy, Jz) is the current 4-vector with p the charge density and J the current 
density and da = = (^^) ^) ^) ^) in (|23] ). We continue to use the Einstein summation 
convention: when an index appears twice, once in an upper (superscript) and once in a lower 
(subscript) position, it implies that we are summing over all of its possible values. So for example 
da *F°'^ implies Yll.=o *F°-^. It is clear that the field Fab is a skew symmetric tensor, namely 
Fab = —Fba, containing six independent real components, the E and B components. We can 
then interpret it as a projective line, since the condition (|52l) is satisfied by plane wave radiating 
fields. In fact, this condition can be expressed as *F°-^Fab = and it corresponds to E • cB = 0. 
Comparing the tensor F"'^ 

I -E^ -Ey 

E^ -cB-, 
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with (|49l) we can write the Pluecker coordinates of the line (BT 



(Ex, Ey, Ez, —cBx, —cBy, —cB^) 



(92) 



which are the E and B cartesian components. 

In spinor form, a real electromagnetic field may be represented by a mixed spinor [33], 

FabA'B' = ^AB ^A'B' + ^AB <^A'B' ■ (93) 

where gab is the spinor metric defined in (l66l) and ^pab is called the electromagnetic spinor. 
Because ga'B' and cab are constant spinors and (pA'B' is the conjugate of ipAB, this explains 
why this (Pab as symmetric spinor encodes all the information. Since cpAB = V^ba is symmetric, 
it has three independent complex components which are related to [33] 



E - jcB 



(94) 
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for real field- vectors E, B. In matricial form, ^pab, a spinor with two indices like eAB is a 
two-by-two matrix. 

The information of the field is bundled in a two-index entity, the electromagnetic spinor ^ab- 
From the quantum physical standpoint, the electromagnetic field is carried by photons, particles 
of spin equal to one. If we want to represent a quantum field with a spinor, a simple rule for 
the number of spinor indices is that the number of indices of like type is twice the quantum 
spin. So ifAB automatically represents a spin 1 boson such as a photon. It is also the case that 
a spinor with primed indices (pA' b' represents the antiparticle field, and in case of a photon the 
opposite helicity. 

Our aim is to attach a physical meaning to the spinor with one index of section |ll] that is to 
derive the polarization state in the form of spinor with one index ip^ that is usually designated 
as the coherent polarization state. However, a spinor with one index, like i/j^ would represent 
quantum fields for fermions with spin equal to |, e.g. the massless neutrino. Therefore from the 
physicist's point of view it may appear manifestly incorrect to represent an electromagnetic field 
by a spinor with one indejq^. 

The solution to this problem, which explains why the polarimetric notation of transverse 
wave states and Stokes vectors has not to our knowledge been formalized in spinor form, is 
that there is missing structure. Absence of this missing structure from the representation means 
that polarimetric tensor and spinor expressions do not appear to transform correctly. From the 
polarimetrist's point of view, absence of the missing structure leaves room for ambiguity in 
the way Jones vectors should be transformed when the direction of propagation is variable, 
something already highlighted in the work of Ludwig [45]. 

In order to obtain the polarization state in one index spinor and see where the missing structure 
is hidden we need to consider the electromagnetic potential as the primary object rather than 
the field tensor. 

B. The electromagnetic potential and the wave spinor 

The expression (|931) is not suitable to derive a one-index spinor to represent a harmonic 
polarization state since the polarization information is equally contained in lpab and ^pa'B'- For 



except in a two dimensional representation [44]. 
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this representation, a more convenient spinor to use, which carries all the required information is 
the electromagnetic potential^. It is a 4— vector whose components are the electrostatic potential 
and the magnetic vector potential A: 

$" = (0,A). (95) 
The electromagnetic tensor Fah can be expressed in terms of $a as [34]: 

f E = V0 

Fah = da^h-dh^a ^ { (96) 

I B = V X A 

where = (0, — A) . If we restrict attention to the Fourier domain then the derivative becomes 
simply an algebraic operation. In fact, 

da ^ jka (97) 
and the (l96l) becomes the skew symmetrized outer product, 

Fah=j{ka^b-h^a). (98) 

Note that, from now on, since we are working in the Fourier domain, the harmonic analytical 
signal representation is implicit. The fields and the derived objects become complex. 

At this point we can notice that, given the potential, in order to derive the field, Fah, it is 
necessary to assume the wave vector, ka. However, in polarimetry, this is the goal: to strip off 
assumed or known quantities: frequency, direction of propagation, even amplitude and phase, 
to arrive at the polarization state. Now, in spinor form, the vector potential is a 4— vector 
isomorphic to an hermitian matrix, ^ab' like any covariant 4— vector (|42l) 

\-A, - jAy <j) + A,j 
Moreover, as is well known, the vector_potential has gauge freedom, so is not uniquely specified 



for a given field. In fact, the field Fai 



is not altered if the potential $a is changed subtracting 



the 4— gradient of some arbitrary function x- 

^a = ^a-daX < ! ■ (100) 

[ A = A + Vx 

^^The vector potential is often used in antenna tlieory, where the potential in the far field can be simply related to each current 
element in the source. 

^'^and consequently E and B. 
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In relativity texts, a particular gauge, namely the Lorenz gauge, is typically singled out, as it 
is a unique choice that is invariant under Lorentz transformations. However, if this generality is 
not required other choices are possible: 

• the radiation gauge, also known as the transverse gauge, has the consequence that 

k^Fat = 0, (101) 

• the Coulomb gauge can be expressed as 

iOa^" = 0, (102) 

with iOa = (1, 0, 0, 0), and implies = 0. 
These two choices are not mutually exclusive and can be simultaneously satisfied for a radiating 
plane wave. Then we have for a wave propagating in the 2— direction that the radiation gauge 
condition implies that = and together with the Coulomb gauge we obtain 

$01' 



AB' 



$10' 



(103) 



-A.^ + J Ay 

-Ax — jAy 

This contains two generally non-vanishing components that transform conjugately with respect 
to one another, and can be identified with components of opposite helicity (the two circular 
polarization components), since the E vector is algebraically proportional to the vector potential 
in the Fourier domain (l96l) . In order to obtain the traditional 'Jones vector' representation as the 
one-index spinor tpA, it is necessary to find a form of projection of ^ab' onto a spinor with one 
index. The polarization information contained in (11031) can be amalgamated into a single spinor 
simply by contracting with a constant spinor 9^': 

■1 ^ nB' ( "^OlA M 

ipA= jko<l>AB'0^ =jko \ = I (104) 

\^$io' J \lj \3ko<l>io'J 

with kQ = - . This achieves the stated goal and jkQ $01' and jk^ $10' are the circular polarization 

componentO- It has been necessary to introduce extra structure, namely the contraction with 6^' 

^'in fact in the Fourier representation, tlie E vector is algebraically proportional to the vector potential A. Using i96l and 
l |97| l we obtain E = —jkoA. Then, jko ^gi/ = — jEy and jko $10' = + jEy are the left and right circular polarization 
components. 
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to obtain the one-index representation tp^. This explains the apparent physical inappropriateness 
of the one-index spinor representation. While the extra structure may for obvious reasons be 
unattractive to the theoretical physicists, it is by contrast of value to the practical polarimetrists 
because it can be inserted when needed to resolve questions relating to the geometry of scattering. 

We have to notice that the choice of the constant spinor 6^' in (|104l) has a degree of 
arbitrariness. In fact a rotation by an angle a about the z— axis is given by (f82l) : 



r = U^;,,e^' = . (105) 

\^e-^t e^'j 

Using the new 6^' in the defining relation (|104l) . the new spinor will be 

[ e-^t ^i'$oiA 

i^A=jko , (106) 

which shows that a rotation around the 2— axis will result in a differential phase offset between 
the two circular polarization components $oi' and $io'. However, if unitary property are to be 
preserved then 

|^o|' + |^i|' = |V'o|' + |^i|' (107) 
which projectively means that the components of 9^' have to be equal in amplitude: 

\e°'\ = \e^'\. (108) 

Spinors with equal amplitude components necessarily correspond to real points 

on the equatorial planeQ. For reasons that will shortly become clear, we call 6^' the phase flag. 

Whilst (11031) employs a special choice of coordinates, it is important to note that a rotation 
of reference frame may be effected through the relation (see (l78l) ) 

^AB' = Ua^Ub?'^cd'. (110) 

The covariant (coordinate independent) formulation of (11041) means that it is valid for waves 
propagating in arbitrary directions, namely it is valid in any rotated frame if all the elements are 
transformed according to the appropriate rule: 

^A = f//V^B, (111) 

^*Note that this fact is valid in circular basis. 
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^AB' = U/Us?'<l>CD' (112) 

§a> = U^^'9b' (113) 

with the unitary spinor describing the rotation of the frame which components can be 
parametrized as in (|80l) . 

This is the most important result of this section: we have been able to define the spinor 
containing the polarization information from the electromagnetic field using a constant spinor 
we named phase flag. This definition is valid in any reference frame, namely it is valid for any 
direction of propagation considered. Table U summaries the main steps to follow to extract the 
polarization information from the electromagnetic tensor. 

V. Spinors, phase and the Poincare sphere 

So far we have deliberately avoided more than scant reference to the Poincare sphere. The 
reason for this is that by integrating the material from the previous section with the projective 
interpretation of Sections HI] and Ull] we are now able to arrive at the first remarkable consequence 
of this work, namely that the Poincare sphere may be identified with the wave sphere by a direct 
geometrical construction. In order to see this, we have simply to interpret geometrically the 
algebraic relations of the previous section. The vector potential in its covariant form may be 
represented projectively as a plane as we have seen in the section [III If no gauge condition is 
specified, then would represent a general plane in projective space. However, the radiation 
gauge condition k^-Fab = (|101l) implies that 

A;"$^ = ^ "^(j) - KA^ - kyAy - k,A, = Q (114) 

namely that is any plane that passes through the point of tangency of the wave plane ka with 
the wave sphere Vt"'^ as shown in Fig. [5] (see (|4T1) and following). The Coulomb gauge condition 
cUa^" = (|102l) . however, imposes the condition that $a passes through the center of the wave 
sphere uJa = (1, 0, 0, 0). Taken together, these conditions imply that the plane $a passes through 
the axis of the sphere normal to the wave plane ka (Fig. (6]). The appropriateness of the projective 
interpretation is seen by noting that the significance of the line intersection of the planes and 
ka is that it represents projectively the electromagnetic field tensor Fab- 
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Let us suppose that the field tensor Fab is given. Algebraically, using homogeneous coordinates, 
Fab (l98] ) in the Fourier domain is proportional to 

Fab (X ka^b - h^a (115) 



and it expresses the dual homogeneous coordinates of a line (I56l)-(l58l). called Pluecker coordi- 
nateo in the projective space [19]. 

Now, given the line Fab and the wavevector ka we have that is a general plane in ¥'^. 
Imposing the radiation gauge condition (|101l) we have that the plane passes through the point 
fc" but it has still the freedom to pivot about the line Fab (see Fig. [5]). The unique representation 
is then obtained by choosing the one element of the family of planes $a that passes through the 
center of the sphere, namely the plane which satisfy the Coulomb condition (11021) . Very nicely 
the geometry neatly illustrates the gauge freedom. 

We now come to our main objective: to express the one index spinor tpA geometrically. In 
the geometric interpretation, the totality of wave states with wave vector A;" are obtained by the 
lines Fab that pass through the k"' axis in the plane ka- Note that we may admit complex lines 
as the geometry is generally valid over complex coordinates. For every line Fab, the plane 
is thus uniquely determined given the gauge fixing. 

Now we consider the spinor equivalent of the plane which is as in (l99l) still 

geometrically represented as a plane. As seen in the section |IIl spinors can be interpreted as 
complex lines generating the wave sphere (see equations (f42)) -(l46l)). We suppose that one of the 
two generators of the wave sphere 6^' is chosen as a fixed reference generator. Then provided 
6^' does not lie in the plane $ab' (which is guaranteed by the gauge choice), there is a unique 
point of the generator 9^' that intersects for any wave state and thus we can write 

1pA = ^AB'0'''. (116) 

The spinor ipA represents the generator of the complementary family to 9^' of the sphere that 
intersect ^ab' in the same point as the —generator, since 

^'in general Fab is a bivector [46] because of its antisymmetry property. For a plane wave radiating field the bivector becomes 
simple and it is representable as the line U15t . The condition for the bivector Fat to be simple is FabF"''' — which is equivalent 
to E ■ cB = and to ( 152b . The condition for the bivector to be null is FatF'^'' = which is equivalent to |Ep = |cBp. This 
condition states that the bivector intersects the sphere 57"''. Both conditions are met by a plane wave radiating field. 
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• a linear relation must exist between these objects, due to the fact that generators of opposite 
type have a unique intersection. In fact, reconsidering the equation (11161) we can rewrite it 
as 

^AB^^O""' = 0, (117) 

having contracted both sides of equation (|116l) with and using i^Aip^ = as in (fTTI) . 
Now, if represents a plane, then 'tp"^9^' represents a point p". This point is in general 
complex since the matrix obtained P^^' = ip^O^' is only singular but not necessarily 
Hermitian as in (l37l) or (|3TI) and it is the intersection of the two generators (two lines) of 
the wave sphere of opposite families, ^/^^ and 6^' . Moreover, this point lies on the wave 
sphere since the matrix P^^' is singular. Finally, the relation (|117l) tells us that that this 
null point belongs to the plane 

• (II 161) is the only such relation that transforms covariantly and homogeneously. 

Now, the reason we claim that we can identify the wave sphere with the Poincare sphere is found 
by considering the effects of pure rotations about the k'^ axis for any wave state: it is seen that 
keeping 9^' fixed when the plane rotates by 180°, the same point of intersection arises. This 
explains the double rotation phenomenon of Stokes vectors, when rotation takes place around 
the wave vector. In fact we can construct the Stokes vector noting that only one point of any 
generator is real. It is clear that 

7/;^^^, = ^^^,= I \ ^ Sa (118) 

is a singular and Hermitian matrix (l37l) representing a coherency matrix for a single wave-state. 
It is isomorphic to a real vector Sa, the Stokes vectoiE! 

S = (,So, 51,52,^3)= (119) 
= (l^oP + |^iP,2Re[^o^i,], -2Im[^o^i,], - l^iH 

whose discriminating condition 

Si - Si - Si -Sl = Q (120) 

^'^The components of the Stokes vector appear to be different from the usual ones, since we have used throughout this paper 
circular polarization basis instead of the usual linear basis as stated in the footnote [3] 
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recognizably represents a purely polarized state. Geometrically the vanishing of the determinant 
of ^AA' expresses the condition that a plane should contain two conjugate generators (^pA and 
^jJA') of the polarization sphere, so that the covariant Stokes vector Sa represents a real tangent 
plane of the Poincare sphere which now is fused with the wave sphere. In fact the condition 
(|120l) can be expressed as 

San''''Sb = (121) 



which is the equation of a sphere in homogeneous coordinates (as in (|44|) and P3l)). 

We finally remark that S is not a true tensor object because, like tpA it omits structure (the 
phase flag essentially) that also has to transform geometrically. 

The introduction of the phase flag is the key element to identify the polarization state for a 
wave vector in any direction with a spinor. Any transformation of geometric orientation does 
not change the spinor algebra. The spinor expression (II 161) transforms covariantly so that the 
relation will remain valid in any rotated frame if all the elements composing the relation are 
transformed with the appropriate rule as shown in (|1 1 ll) - (ll 131) . The essential distinction between 
geometric and polarimetric basis transformations is that in the former case the phase flag 6^' 
must be transformed for consistency, while in the latter 9^' must be regarded as fixed. Fixing 
9^' and varying ^ab' we obtain all the polarization states for one direction of propagation. 

VI. Illustrative examples 

We will now report some numerical examples in order to illustrate how all this really works. 
We start first to choose one direction of propagation, which is as usual the 2;— direction. The 
covariant 4— vector ka is then 

A;„ = A;o(l, 0,0,-1) (122) 

with ko = ^. We choose now a linear polarized wave at 60° to the x— axis propagating along z. 
The electromagnetic field tensor will be 

/ cosf sinf \ 



ab 



COS I COS I 

- sin I sin | 



(123) 



\ -cosf -sinf / 
where E is the amplitude of the wave, cut — k ■ x = is the phase (|22l) and a is an initial phase. 
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In order to isolate $a from Fab, we use the (|981) contracting Fab with = (1,0,0,0). We 
obtain 

Fabt' =j{ka<^b-kb<^a)t'' = -J $a (124) 



and has components 



A;o V 3 3 

Using (|99l) . the corresponding is 



$^ = ^ ^e^('^+«) ("o, - cos ^, - sin ^, ) . (125) 



?■ V , \ / — cos ? + 7 sin ^ 

= _L^e^('^+") 3 ^ I . (126) 

^0 \_cosf-jsinf 



Now after selecting a convenient fixed generator 9^ 

V2\l 



(127) 



we are ready to compute the polarization spinor ipA with the (|104l) : 

i^A = jko^AE'e""' =jko{^AO'0'' + ^Al'e'')= (128) 

, 1 / e^-'t 

which represents the Jones vector in circular polarization basis for a 60° linear polarization. The 
corresponding polarization ratio is: 

^ = e^'^. (129) 



We can calculate the corresponding Stokes vector using (II 191) and we obtain: 

/ 27r 27r \ 
5„ = £;Ml,cos— ,sin— ,0 . (130) 



3 3 ^ 

which is the Stokes vector for a linear polarization with orientation angle tp = 60°. 

The next step is to see what changes if we choose another phase-flag 9^' , keeping in mind 
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we obtain 

V'A = jko^An'O""' =jko{^AO'0'' + <^Ai'e'')= (132) 

The corresponding polarization ratio and Stokes vector are 

fi = e2^(t-f ) (133) 




Sa = [I, cos2(- - g ), sin2(- - g), Oj • (134) 

which is a linear polarization with orientation angle ^ = 15°. It is clear that the phase flag 
9^' keeps trace of the reference for the orientation angle of the polarization ellipse (see Fig. 
HI). We can notice that the components of the phase flag 9^' apply phase factors to the circular 
polarization components as we have shown in (11061) . We can conclude that this phase offset is 
relative to the origin of the orientation angle. The phase flag is the element which defines the 
orientation angle and the phases of the circular polarization in a consistent manner. For these 
reasons it is the key to define polarization states for arbitrarily oriented wave vectors. 

Now we can do the same calculation using the tensors instead of the spinors and for brevity we 
omit the multiplication constants. We have the plane $„. We calculate the Pluecker coordinates 
9°-^ of the line corresponding to the generator 6^' and then the intersection of the plane potential 
with the line generator, in order to find the polarization state, the generator of the other type. 



We can use (1461) where A^' = ( , ) and /€a = ) • We obtain: 





^^a\b' =11 = 1^'' ■ • (135) 

-X - jy 

Since we want to compute the generator through 9^' then using (11271) A = — 1. The corresponding 
one index tensor is: 

r = (136) 
In order to compute the projective Pluecker coordinates of this line, we consider two points on 
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this line for example = (1, 1, 0, 0) and q"- = (0, 0, j, —1), and we use the relations (|49l)-(l5T 

/ j -l\ 

J -1 

-J -J 

\ 1 10 0/ 

(0,-j,l,0,l,j). (138) 



6' 



(137) 



In order to find the generator of the other type, the polarization state, we compute the intersection 
of the line O'^^ with the plane $b: 

The spinor corresponding to this point will be computed using (|135l) : 

T + Z 



(139) 



(140) 



with of course r = sin |, x = sin |, y = — cos |, z = —j cos |. The spinor can be expressed 
as: 

^ 1 



e' 3 



(141) 



which is projectively the same as (|128l) . 

Now we apply a rotation in space to the bivector F"''', for example a rotation of 45° about 
the ?/— axis. The rotation matrix turns out to be: 

/l \ 



cosf 




\0 sin f 

The new electromagnetic field tensor will be: 

Pcd 

( 



1 





— sm 


COST 



(142) 



pC pa 77IC1 



(143) 



cos J cos 3 



sm 3 



cos J cos 3 



sm 



sm J cos 3 







— sm J sm 3 



cos 3 



sm J sm 3 



— cos J sm 3 



sm J cos 3 



cos 3 



cos J sm 3 
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We can notice that the equivalent matrix calculation to (11431) would be: 

F = RFR^ 



(144) 



where F, F and R indicate the corresponding matrices and R^ is the transpose of the rotation 
matrix R. 

The new wave vector will be: 

'd _ r>d la 



and 



A;'^ = i?'*^r = (l,-sin^,0,cos^ 



TT TT 

1, sin — , 0, — cos — 
4 4 



The new potential can be computed like in (|124|) and we obtain: 



7r TT TT TT TT 

= 0,-7 COS — cos — , — 7 sm — , — ? sm — cos — 



(145) 



(146) 



(147) 



We can easily verify that the potential is still a plane through the origin uj"- and through the 
point k"--. 

^aOJ"" = 0, 4> J" = (148) 

In this calculation we are changing the orientation and calculating the corresponding new 
polarization state. In order to compute the new polarization state we have to rotate the generator 
6'^'^ as well, to obtain: 

9"'^ = R^R^e""^ = (149) 

/ sinf j -cosf\ 

-sinf jcosf -1 
-3 -jcosf -jsinf 
\ cos f 1 3 sin f / 

In order to find the polarization state we need to find the point where the line 6^'^ intersects 
the plane $a as before^!. 



_ (jab^^ _ ^gj^ _^ ^ _ _ _|_ _ ^^g _ 

TT TT TT TT TT 

— cos — , — 7 cos — cos — h sm — sm — 
3 3 4 3 4 



TC 



TV TT 



TV TV 

— cos — . 

3 4' 



(150) 



Of course = R'^^^'i but we have gone through the full calculation again in order to make clear how tensors and spinors 
work. 
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Using the definition (11351) we can calculate the polarization ratios of the two generators: 



^_ ^ + ^ cosg -smg ^ ^^^^^ 

-X + jy cos I + e^if sin | ' 

A= " + ^ = cosf + sinf ^^^^^ 
—X — jy — cos I + sin I 

Now we perform the same calculation using spinors. We want to calculate the new polarization 

state for a new orientation k"-. It will be much faster and simpler. The unitary spinor corresponding 



to R% is Uj^: 



TT A A ■ • / A A \ 

U B = ~ cos — e ^ — J sm — (r ujb + uj tb) 



cos ^ sm 



sm I cos g 



(153) 



where 




are the spinor corresponding to the axis of rotation, the y— axis Ta = (1,0,1,0) and the 
orthogonal one Oa = (1,0,-1,0). represents the rotation in space the generates the new 
orientation A;". 

The new polarization spinor can be simply computed as 

~A 4 R je~^^ / — e-'^ cos f + sin f , 

= Uj,r = ^ I _ L _'\- (155) 

The corresponding covariant form of the spinor is 



cos I + e-^ 3 sm ^ 



? e"-^ 3 / — cos 5 — e-'^s'' sm ^ . 
i^A = ^-^ .J _ ' (156) 



2w 

3 cos ^ + sm ^ 



V2 

and the corresponding polarization ratio is exactly (11511) 



1— TT - TV 

& 3 cos -J- — sm o 
cos I + e-' 3 sm I 



The new phase flag 9^' can be simply computed as 



1 / cos f + sin 5 

7! 



r = [/ir = ^ r « ' ' I . (158) 



cos g — sm - 



41 



Using tlie conjugate version of (|68I) . the corresponding covariant phase flag is 

= 4/ 1 / - cos f + sin f \ 
Ob' = eA'B'd^ = ^ ' M (159) 

V 2 \ + cos f + sin f / 



and the corresponding polarization ratio is (11521) . 
The corresponding Stokes vectors are 



•0 ^ = ( l,cos — cos-,sm— ,cos — sm- ) , (160) 



e^' Gb = (^l,-cos-,0,-sin-j . (161) 
tp"^ is the new polarization state for a new orientation k"- which corresponds to a new phase flag 

QB' 



VII. Discussion and conclusions 

The representation that has been arrived at allows a spinor to describe a polarization state 
for a wave vector in any direction. In any fixed direction, we extend to Jones vector calculus 
in any chosen basis by applying unitary transformation to the polarization spinor alone, and 
not to the phase flag. The natural representation turns out, as might be expected, in terms of a 
circular polarization basis. If we change the direction of the wave vector we have also to apply 
the unitary transformation to the phase flag. 

Geometrically this representation can be determined since two generators of the sphere iJja 
and 6^' (one of each kind) pass through any point of it, and then through every tangential plane 
representing a wave vector there pass two generators tpA and iJja' which are complex conjugates. 

For the coherent field, one represents the polarization state of the propagating wave as 

i/jA a e^'^-^^" = e-'^'^*'*' '') (162) 
and the other, the conjugate field as 

ijA' cc e^^'^"^" = e^(-^i+^-^), (163) 

These fields are conjugate solutions, but both propagate in the same direction, as the equations 
for constant phase surface are identical. The use of strict spinor algebra prevents the often- 
committed error of associating a conjugated Jones vector with a backward propagating field. 
In each wave plane, therefore, there is one generator for the unconjugated Fourier component. 
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and it is possible to show that the entire collection of these forms a ruled surface (regulus) 
on the wave sphere. The generator lying in the plane {ko, k) corresponds to the wave state of 
— 1 helicity (LHC), while the generator in the plane corresponding to the 'backward' direction 
(fco, -k) is that of +1 helicity (RHC). 

This is the origin of the apparent conjugation/reversal symmetry. The lines on the regulus 
form a one-dimensional linear space, and in this sense we can say 



for a homogeneous projective spinor in which the complex polarization ratio ^ = ^ runs from 
to oo. As we have shown there is an isomorphism between the generators of the sphere and the set 
of the spinors ^pA, and this can be effected in an invariant way with respect to any linear change 
of basis. Thus we have arrived at a unified polarimetric description in which both the geometry of 
'world space' and the abstract mapping of the generators of the sphere are handled consistently 
using spinors. The geometric interpretation we have introduced explains the fundamental place 
of the Poincare sphere in polarimetry as an invariant object under linear transformations with 
its reguli whose generators are wave spinors, constituting its invariant subspaces. Identifying 
generators as states of polarization, the structure of the Poincare sphere is preserved under all 
linear processes. The sphere is considered an invariant of the theory, the 'absolute quadric' [19]. 
The Poincare sphere and the wave sphere are hereby unified. The well-known phenomenon 
of 'double rotation' of Stokes vectors with respect to rotation of world coordinates is down 
to the fact that 6^' is not rotated when a basis transformation is made while for geometric 
transformations the phase flag is included (|113l) . It may appear that the phase flag, 9^' was 
introduced in a somewhat ad hoc manner, and indeed the choice 



was not a unique one. It is obvious that the absolute phase of 9^' could be chosen arbitrarily and 
even that its elements undergo relative phase transformation without loss of information. This 
is indeed the case, and so it appears wholly reasonable that we name this object the phase flag. 
Its practical importance is that it defines a phase reference for both components of the polarized 




(164) 




(165) 
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wave. This is to be contrasted with conventions fixing one component as real, eg 




(166) 



which fails when 9 |. In a further planned paper on this topic we shall address the problems 

of defining the phase of an electromagnetic plane wave, and identifying the role of the phase 
flag in the relationship between geometry and phase. 
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TABLE I 



The electromagnetic field 



Fab = jika^b - h^a) 



ab 








Ex 


Ey 


\ 




-Ex 








CBy 




-Ey 








-cBx 


V 





-CBy 


cBx 


/ 




projcctivcly a line 
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The electromagnetic potential 
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The polarization spinor 
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R P Q 
• • 

Fig. 2. A "real" visualization of the complex line. P = R-\- (^Q is any point on the line. For = ((^ = 0), the point P 
coincides with R, and for — {(^ — oo) P = Q. 
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Fig. 3. A quadric surface in the form of a hyperboloid. On the left we can see the hyperboloid generated by one line rotating 
around one axis. In the picture on the right we show the two family of generators. 



y 




Fig. 4. An example where the same linear polarization is defined using two different phase flags. In one case we have 
used 9^ = (^) which is the x— axes in circular basis. The orientation angle is ip = 60°. In the other case we have used 
6^ = [ e+^?) w^'i'-h is ^ °f respect to the a;— axis. The resulting orientation angle is -i/i = 15°. 




Fig. 5. Relation between wave vector, electromagnetic field and the potential for the radiation gauge. In this gauge, the plane 
$a is any plane intersecting the point fc" on the wave sphere $1"'', hinging around the line Fab which is the intersection between 
the plane and the wave plane ka- 
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Fig. 6. Relation between wave vector, electromagnetic field and the potential for the Coulomb gauge. In this gauge, $a is the 
plane through the point fc" and the center of the wave sphere £7"''. The electromagnetic field is the same line of Fig. |5] the 
intersection between the plane $a and the wave plane ka. 



